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Abstract 

Water and ice Cherenkov telescopes of the present and future aim for the detection of 
a neutrino signal from extraterrestrial sources at energies E u >PeV [Woschnagg and 
AMANDA Collaboration (2004); Montaruli (2003); IceCube Collaboration (2004)]. 
Some of the most promising extragalactic sources are Active Galactic Nuclei (AGN). 
In this paper, the neutrino flux from two kinds of AGN sources will be estimated 
assuming p7 interactions in the jets of the AGN. The first analyzed sample contains 
FR-II radio galaxies while the second AGN type examined are blazars. The result 
is highly dependent on the proton's index of the energy spectrum. To normalize 
the spectrum, the connection between neutrino and disk luminosity will be used 
by applying the jet-disk symbiosis model from Falcke and Biermann (1995). The 
maximum proton energy and thus, also the maximum neutrino energy of the source 
is connected to its disk luminosity, which was shown by Lovelace (1976) and was 
confirmed by Falcke et al. (1995). 
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1 Introduction 



Large v olume neutrino te l escopes like AMANDA, Antar es and IceCube |Woschnagg and AMANDA Col 
( 153) ; iMontarnlil gfloj ; iTceCube Collaboration! jjjpl] are being built to de- 
tect the neutrino flux from extraterrestrial sources. For a diffuse analysis, the 
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detection of such a signal is possible at energies at which the extragalactic neu- 
trino flux contribution exceeds the number of atmospheric neutrinos. Active 
Galactic Nuclei (AGN) are believed to make up a lar ge fraction of the extra- 
galact ic neutrino flux. While previous models - see e.g. (Learned and Mannheim 
( 200(1 )] for a review - used the Cosmic Ray flux or a certain part of the photon 
spectrum to estimate the neutrino flux from A GN, the model presented in 
this paper uses the jet-disk symbiosis model by Falcke and Biermann ( 19951 ) 
and allows for the dependence of the maximum particle energy on the source 
properties. This is the key point of our analysis. 



The total neutrino spectrum at Earth is given as 

dN , „ n , N dn , r s dV 

dz dL 

z L 
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The parameters in the calculation are the integral radio luminosity L and the 
redshift z. The energy of the neutrinos at the source E v is redshifted due to 
adiabatic energy losses on the particles' way toward Earth, E v = (1 + z) ■ E®, 
where E® is the neutrinos' energy at Earth. dN/dE v is the single source spec- 
trum of an AGN, while dn/dL is the radio luminosity function (RLF) per 
comoving volume dV/dz and luminosity interval. To obtain the total number 
of AGN per redshift interval, the RLF has to be weighted by the comoving 
volume and divided by a factor of knd 2 L with dL as the luminosity distance to 
receive the neutrino flux at Earth. Finally, the result has to be integrated over 
all possible redshifts and luminosities, where the integration limits depend on 
the source sample as discussed in section 4. 

The neutrino flux will be calculated for two different source samples, one in- 
cluding FR-II sources with strong extended radio emission (in the following 
referred to as "steep spectrum sources") and the second one flat spectrum 
blazar sources (referred to as "flat spectrum sources" throughout the paper). 
It will be assumed that neutrinos are produced by pions resulting from a A- 
resonance of the initial p 7 interaction. In this scenario, the neutrino spectrum 
follows the proton spectrum. Multipion events (see e.g. 
will not be included in this calculation, since the changes of the spectrum 
due to these events would lie within the estimated uncertainty of the calcu- 
lation due to uncertainties in the jet-disk symbiosis model and the sources' 
Radio Luminosity Functions (RLFs), see section 5.2. For the steep spectrum 
sources, the spectral index of the proton spectrum will be calculated from the 
synchrotron emission of the electrons. This is possible since the proton and 
electron spectra resulting from shock acceleration are equivalent. This method 
cannot be used for the flat spectrum source sample, since the radio emission 
of the jet is spatially unresolved and therefore, the observed spectrum is a 
superposition of the different radio spectra. In this case, the spectrum will 
be a superposition of the spectra of all knots, each with a particle index of 
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2, following the d erivation of the index in diffuse shock acceleration, see 



Protheroel ( 19981 )]. The maximum energy of the proto ns in the j ets w ill 



p 2 
e.g 

be assumed to be luminosity dependent as is suggested in Lovelace! ( 1976b ]. 
To norm alize the generic AG N energy spectrum, the jet-disk symbiosis model 
given m |Falcke et all (Il995h ] will be applied to connect the neutrino luminos- 
ity with the disk luminosity and therefore, also with the radio luminosity of the 
jets. Throughout this paper, an Einstein-de Sitter Universe (Q m = 1, Q\ = 0) 
with h = 0.5 is assumed, since the used RLFs are not available in the con- 
cordance mode l (fl m ~ 0.3, fl\ ~ 0.7) which has been confirmed by recent 



measurements [Spergel et al. ( 20031 )] . However, the resu lt should be indepen- 
dent o f the cosmology used as it has been pointed out in |Dunlop and Peacock 

The paper is organized as follows: In section 2, the two AGN samples will be 
discussed. They will be used to determine the AGN radio luminosity function 
(RLF). In case of the steep spectrum sources, the spectral indices between 2.7 
and 5 GHz are used to evaluate the index of the primary particle spectrum. 
Section 3 gives an explicit description of the normalization of the single source 
spectrum to the disk luminosity of the objects. In section 4, the integration 
limits for the redshift and luminosity distribution are discussed. Finally, the 
different models for both types of AGN sources are applied to calculate the 
AGN neutrino flux. It will be shown that the result strongly depends on the 
index of the particle spectrum. 



2 The samples 



The fi rst sample consists of 356 steep spectrum s ources, selected bvL Willot t et al 
( 200 1|) . The sourc es are s ampled from the 7CRS McGilchris and Rilevl |l990|) : 
Lacv et all{l99l]. 6CE (Rawlings et al.l |200ll)] and 3CRR |Laing et al.ldl983h] 
catalogs with flux measureme nts at 178 MHz. The second sample comprises 
171 flat spectrum sources from Dunlop and Peacock ( 199fj ). They are selecte d 
from the Parkes sele cted regions sample, the Parkes ±4° zone IPeacock ( 1985 ) : 



Wright et all (|l982h ]. the ' Northern-Sky 1 suryev of IPeacock and Walll ([l98lT ) 



and the 'All-Sky' survey of IWall and Peacock! (|l985T ). 



2.1 The spectral and particle index of the sources 



As discussed in detail in, e.g., |Rvbicki and Lightman ( 1979! )]. the index p of 
the proton spectrum is correlated to the index a of synchrotron radiation by 
p = 2a + 1. 
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2.1.1 Steep spectrum sources 



23 of t he total 356 sources are identified in the S1-S5 catalogs | Pa u inv-Tot h et al 



Paulinv-Toth and Kellermannl ()l 972h : IPanliny-Tot,h et alJ (|l 978h : lKfihr et al 
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( 198ll )] which give the spectral index of the sources at 2.7 — 5 GHz. To estimate 
the steep spectrum index a s at 2.7 — 5 GHz, the maximum in a distribution 
of the source counts versus a s is determined applying a Gaulfi fit (see Fig. 1). 
The maximum is found to be at a s = 0.8 with a peak width of o~ s = 0.2. In 
calculations of the AGN evolution functions by Willott et al., a spectral index 
of a s = 0.8 is assumed | Willott et al.l (|200l| )]. The spot check of the sources 
supports this a ssumption. Hence, the index of a = 0.8 is adopted. This result 
is supported bv lGregorini et all (|l984l ) who analyze the spectral index of radio 
galaxies in the range of 2.7 to 5 GHz resulting in an average index of around 
0.8 for steep spectra. Consequently the particle index is p s = 2 • a s + 1 = 2.6. 
Here, it is assumed that the stationary particle spectrum is equal to the par- 
ticle injection spectrum, which is based on the assumption that particles from 
this region undergo energy scale independent losses such as would occur in a 
flow. There is unfortunately no good and well-tested model available at this 
time for the particle transport, convection versus diffusion and other processes, 
so that we make this simple assumption by Occam's razor. To demonstrate 
the influence of the particle index on the resulting spectrum, the calculation 
will also be done with different harder spectral indices. 




a. 



Fig. 1. Histogram of the spectral indices of the steep spectrum sources. Sources 
counts are shown against the spectral index at 5 GHz. The maximum is found at 
a s = 0.8. The width of the Gaufi distribution is a s = 0.2. 
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2.1.2 Flat spectrum sources 

Since the radio spectrum of blazars is spatially unresolved, the observed in- 
dex of the radio spectrum cannot be identified with the index of synchrotron 
radiation. The particle spectrum is a superposition of the spectra of the single 
knots, 

dN/dE v ocJ2 a i E ~ P 



assuming the same particle index for all sources, p ~ 2, as suggested by diffuse 
shock acceleration for strong shocks, and also by the radio data of strong 
jets. An index clo se to y ~ 2 is supported by the individual radio spectra as 



described in, e.g., (Bridle and Perlevl (|l984)] for strong sources. 



2.2 Radio Luminosity Functions 



2.2.1 Steep spectrum sources 

The AGN RLF depends on the luminosity and on the redshift. To find the 
proper relation empirically, a factorizing separation of the density in a lumi- 
no sity dependent funct ion and a redshift dependent distribution is assumed 



bv IWillott et alJ (j200lh . The model includes the steep spectrum sources as 
they are explained above. The AGN RLF is given at a frequency of 0.151 GHz 
depending on the differential luminosity Lo.isi- The RLF was assumed to con- 
sist of two separate distributions, a low luminosity function including objects 
without or with only weak emission lines and a high luminosity function with 
objects with strong emission lines. The density p(L .i5i, z) is given as the prod- 
uct of a pure luminosity function pz(Lo.i5i) m units of 1/ (Gpc 3 A log L0.151) and 
the dimensionless evolution function f(z): 

p(L .i5i, z) = p(L .i5i) • f(z) ■ 



Two AGN populations contribute to the RLF: 

Low luminosity sources only show weak or no emission lines. These sources 
are FR-I and weak emission line FR-II galaxies. Above a luminosity L 0A51 ps 
10 25,5 W/(Hz sr), the sources can almost exclusively be assumed to be FR-II 
galaxies. For the low luminosity function, pi(L), an ansatz of 

MU51 



pi(L . 



L 



151, 



Pi 



0.151 



exp 
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is made. This an alytic ansatz is known as the Schechter function in literature 
|Schechte"rl (|l976| )]. a\ is the power law shape, p° is the normalization constant 
and L\ is the break luminosity. The evolution function of the low luminosity 
population is taken to be 



'1 + z) k ' for z < z? 



zf) h for z > zf 



The redshift evolution of the sources is known up to ~ zf. At higher red- 
shifts, it is assumed to be constant since there are experimental indications 
Mivaji et al. 1 200(1 )] of a constant or slightly decreasing distribution at red- 
shifts up to 2 ~ 6. The values of five free parameters used in the functions 
are listed in table 1. The complete (luminosity and redshift dependent) RLF 
is shown in Fig. 2. 

The high luminosity population consists of FR-II galaxies with strong emis- 
sion lines. The form of the high luminosity function, Ph(L i 51 ), is similar to 
the one for the low luminosity population. Here, however, the exponential 
function applies at low luminosities while the power law dominates at higher 
luminosities: 



Ph(L 



0.151, 



P°H 



^0.151 



-Cih 



■ exp 



L 



0.151 



The evolution function of the high luminosity population is assumed to be 
an exponential function up to a certain redshift z\ and then continuing as a 
constant: 



fh(z) 



exp 



for z < z\ 



for z > z\ 



A comparison to X-ray data Mivaji et al. ( 2000| )] justifies an exponentially 
increasing evolution. 



The luminosity, which is here given at 0.151 GHz, is used in the integral form 
for further calculations. Thus the luminosity L0.151 is converted to the integral 
radio luminosity in units of 10 42 erg/s, L 42 , by assuming a spectral index of 
a = 0.8. Furthermore, the RLF is given in units of Gpc -3 (A log(L)) -1 . In 
following calculations, the RLF will be given per luminosity interval <iL 42 . 
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Table 1 

Parameters used in the calculation of the AGN RLF for steep spectrum sources 
Willott et ahl fenmh ]. p*} and p° h are given in units Gpc _3 A(log L). L* h and L\ are 
given in W/(Hz sr). 
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log(U) 




25 23.5 24 24.5 25 25.5 26 26.5 27 27.5 2 

!og(L 17 /W/Hz) 



Fi g. 2. Steep spectrum RLF according 
to IWillott et alJ l|200lh versus luminos- 
ity for z = 0, 1, 2. 



Fig. 3. RLF for flat spectrum sources in 
a pure lu minosity evolution model ac - 
cording to IPunlop and Peacockl (jl990h . 
Four different redshift values are used, 
z = 0, 2, 3, 4. 



Flat spectrum sources 



171 fla t sources from four catalogs at 2.7 GHz are used bv lDunlop and Peacock 
( 199(1 ) to determine the RLF of flat spectrum blazars as described above. The 
ansatz for the flat spectrum RLF is a pure luminosity evolution function of 
the form 



p(L,z) 




The function is labeled pure luminosity evolution since only the break lumi- 
nosity L c evolves with the redshift. Here, a and f3 are the power law slopes, 
Po is the normalization and L c (z) is the break luminosity which evolves with 
redshift: 

log[L c (z)] = a + ai ■ z + a 2 ■ z 2 . 
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The luminosity L is given in units of W/(Hz sr) at a frequency of 2.7 GHz. 
The RLF is shown in Fig. 3 with the set of parameters as given in table 2. 
For further calculations, the luminosity at 2.7 GHz is converted to an integral 
luminosity by assuming a flat spectrum with a/ = 0. The RLF is given in 
units of Gpc -3 • (AlogL) -1 , and thus the differential RLF is 




P0 


a 


P 


ao 


ai 


a2 


10 - 85 Gpc" 3 


0.83 


1.96 


24.89 


1.18 


-0.28 



Table 2 

Parameters for the flat spectrum RLF Peacock! 



3 The single source spectrum 



The generic AGN neutrino flux is adopted to follow the proton spectrum from 
the A-decay. This can be described as a power law with an index p having 
an exponential cutoff. This cutoff depends on the particle's maximum energy. 
Also, the power of the spectrum changes by —2 when synchrotron losses of 
pions start to dominate at energies of E h J eak « E™*/&.7. The break energy 
is determined by the ratio of the pion's to the protons mass, m n /m p = 1/6.7, 
since the maximum energy due to s ynchrotron losses of a charg e d par ticle is 
proportional to the particles mass |Biermann and Strittmattei dl987t )]. The 
spectrum is then given as 



dN 



disk ) 



E~ p exp 



Pmax 
J-'v 

break 



for E v < E 



break 



(L disk )-E~p[E u /E b J eak ) exp(--lfe) for E v > E 



nbreak 



.(3) 



The power law behavior is due to diffuse shock acceleration of the protons in 
the jet's shock fronts. The spectrum is limited by the strength o f the magnet ic 
field of the source due to the spatial limit of the Larmor motion Hillasl ( 1984 )]. 



This determines the maximum energy of the accelerated protons, there is an 
exponential cutoff in the spectrum. The magnetic field is connected to the disk 
luminosity of the AGN and subsequently the maximum proton energy can be 
expressed in terms of the disk luminosity, see (Lovelace! (jl97fi! )]: 



E™ x oc 



■'disk 
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as it h as been confirmed by the jet-disk symbiosis model, see iFalcke and Biermann 



BiermannT i 2001 ) , a satura 



( 199fih ]. At high luminosities of L disk > 10 46 erg/s 
tion of the proton maximum energy is expected due to photo-pion production 
in the disk radiation field. This effect will be discussed in section 4. 



The normalization of the generic neutrino spectrum, O , can be found by 
assuming that particles from an AGN produce the same order of magnitude 
of lumi nosity as the j et po wer, Qj e t, and neutrinos contribute with a fraction 



q v < 1 Falcke et al. ( 199R| )]. The observations are best fitted by this concept 



and this number constraint. It is assumed that the power of the neutrinos is 
directly connected to the power of the jet. This in turn is a fraction q > 1/3 of 
the disk luminosity. Here we use the lower limit, q = 1/3, as an estimate. The 
factor q v results from assuming a kinematic factor of ~ 1/2 which accounts 
for the fact that only half of the pro tons are converted i nto pions - about 50% 
of the protons convert into neutrons Miicke et al. ( 2000l )]. Additionally, a (non 
thermal) efficiency of e nt h ~ 10% is applied - 20% of the proton energy goes 
into secondary mesons and 50% of these 20% go into neutrino production. 
The optical depth of AGN sources r p7 can be described as the product of the 
photohadronic cross section <t P7 and the photon density n 7 , and thus 



pi 



^46 



Here, r\i := r/10 17 cm is the extension of the source in units of 10 17 cm. Scaling 
this optical depth to the Eddington luminosity and to the Schwarzschild radius 
demonstrates that the optical depth is independent of the power of the source 
at equivalent distances from the central engine. It also shows that large optical 
depths are easily achieved. The neutrino measurements can then be used to 
set limits on the optical depths, and so on the structure of the emitting region, 
as we show below. Taking into account interactions of protons in Bethe-Heitler 
pair production, an effective opacity for photopion production is reduced to 

T eff = r P1 -V- 



r] gives the reduction of the optical depth due to the optical depth of the 
source taking into account Bethe Heitler interactions, r BH . In all following 
calculations, r e ff = 1 will be assumed. Detailed geometric models for emis- 
sion from zones where the effective optical depth is large, where it is close to 
unity, and where it is much less than unity, need to be considered in a later 
stage of such calculations. This might be easier in the future once there are 
good constraints, or even exact spectra from multi-GeV to multi-TeV obser- 
vations of AGN. Since the neutrino flux is proportional to T e ff, a lower limit 
as an estimate of a maximal optical depth is possible when comparing current 
neutrino flux limits to the calculated flux as it is described in section 5.4. 
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The fraction q u is then 

£ ~ 1 

Qv — q ■ T e ff ■ i^kin ■ tnth ~ — • 

60 



Thus the neutrino luminosity from the jets is given as a fraction of the disk 
luminosity: 



rmin I?min ^ V ' 

^max 



J ^ ■ E^- 1 ■ E b J eak2 exp dE v = q u ■ L dlsk 



J^break 



The integral can be evaluated by developing the exponential function in a 
Taylor series and interchanging the integral with the sum. The lower integra- 
tion limit is a fraction of 1/4 of the rest mass of th e pion , E™ m = E® ± /A = 
(139.57018 ± 0.00035)/4 MeV, (Particle Data OrourJ (|2004 ]. The upper inte- 



gration limit is assumed to be varying with the disk luminosity, see Equ. (8). 
The normalization is thus given as 



max 



x L dlsk = In j + {f(p - 2, 1/6.7, 0) 

+ [f(p, 1,0)- f(p, 6.7, 0)] /6.7 2 } • E ™" +2 

- f(p - 2, Er x /Ef n , -p + 2) ■ Er~ p+2 



with 



oo ( — l) n 

f{a,b,c)= £ , V b -n +a+c 
n -l n\ ■ (n — a) 



The series' are evaluated numerically for each case. 
3.1 Jet-disk symbiosis 



To convert the disk luminosity into the radio luminosity of the jets, the jet-disk 
symbiosis model by Falcke et al. ( 19951 ) is used. 
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3.1.1 Flat spectrum sources 



Assuming a power l aw slope of a 
erg/s, is given as 



Falrke and Biermann 

/ \0.83 
X e) 



the lumin osity of the jet in units of 



U = 6.7-10 42 ^- J D 217 • sin ai7 (z obs 



1.8 



1.42 T 1-42-5 



(4) 
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Parameters appearing in Equ. 4 are given in the parameter list below. In 
using here a = instead of the sum of components with a = 1/2 we follow 
the earlier work. It is not known what the energetic particle spectrum does at 
low energies in the jet. For the final conclusion below we do use p — 2, which 
would correspond to a = 1/2 if pure synchrotron radiation was observed. 



The parameters from above are defined in the parameter list below. The cor- 
relation between disk luminosity and radio luminosity in units of 10 42 erg/s, 



J 42; 



IS 



■'disk 



(2.1 ± 1.9) • 10 45 (L 



\0.79 



12. 



erg 



(5) 



The error has been estimated by taking into account a scattering of the data 
by about a factor three around L{ 2 {Ldisk)- 



3.1.2 Steep spectrum sources 

The model for extended, optically thin sources will additionally be used for 
the calculation of the steep spectrum source flux. The relation between the 
differential radio luminosity at 5 GHz of the lobes and the disk luminosity is 

P -18 10* ^ ( GHZ T ^ /4 -^-<ioo-(^) 3/2 ((]) 
P^e-1-8-10 - m \ c —) p . (6) 

This resu lt is based on the in v estiga tion of a quasar sample, described in 
detail in |Falcke and Biermann ( 199R[ )]. The disk luminosity is connected to 



the integral radio luminosity, L s , as 

L disk = (21.6 ± 16.6) • 10 45 • (1 + z)- 1/2 (Ll 2 ) 2 ^ , (7) 



assuming a spectral index of a = 0.8. Again, the error has been estimated 
considering the scattering of the data by a factor three around the predicted 
result. The parameters which occur in Equ. (4) and (6) are given below. 
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PARAMETER LIST 



Pj-. Jet velocity in units of c which can be assumed to have a value of f3j « 0.986. 
jj: jet's Lorentz Boost factor, jj = 6. 

P-12- Pressure parameter. The external radiation pressure at the AGN is given 
by Pext = P-12 ■ 10~ 12 erg/cm 3 . The pressure evolves with the time due to the 
expansion of the Universe. The pressure parameter depends on the redshift as 
P_i2 = (1 + zf ■ P° 12 , where P° 12 is the pressure of the Universe at z — 0. P£^2 
is set to P®i2 = 1. 

x' el00 : Modified electron density. Let the ratio between the relativistic electron 
density and the total number density of protons be x e and the minimum Lorentz 
factor of the relativistic electron population divided by 100, 7 e ,ioo- The modified 
electron density is then x' elQ Q = 7f7n nZ P . V = 2 is a parameter. x ' e 100 1 has 



been used in all following calculations Falcke and BiermannI (|199Sh 



Qj/l- Total jet power Qj et in units of the disk luminosity, with a given value of 
Qj/i = Qjet/Ldisk = 0.15to;i- 

L^q: Disk luminosity per 10 erg/s: L4Q = L^ s fc/(10 46 erg/s). 

7^5: Jet's Lorentz boost factor, divided by 5. It is assumed to be 7^5 = (6 ± 2)/5. 

U3: Ratio between the total energy density in the je t and the magnetic energy 



density, divided by a factor three. It is set to U3 = 1 [Falcke et all (|199, t Tl )]. Note, 
again, that this fits the observations best. 

£ = 0.15: Fit parameter resulting from the ansatz that the bulk proper velocity de- 
pends on the accretion rate and thus on the disk luminosity, see Falcke and BiermannI 

ET 



iobs = I/73 = 1/6: angle between observer and jet axis. The angle is taken to be 
l/7j, because this is the maximal boosting range and has a solid angle in which 
optimal boosting occurs. We use for simplicity the same angle for all sources: In a 
more refined model obviously isotropy should be used, but then mediated by the 
selection effects, which strongly favor small angles, as are used here. 

D = ^- — cos iobs)' 1 ~ 6: Doppler factor. In this calculation, the Doppler fac- 
tor is assumed to be the same for all sources. A variation of the Doppl er factor with 
the so urce luminosity has been discussed - for an example see e.g. (Falcke et al 



ljl99,^ ]. However, since there is no well-founded model to apply individual Doppler 
factors so that a constant factor has to be used as an approximation. 



Finally, the jet-disk symbiosis can be used to express the generic neutrino 
energy spectrum in terms of the radio luminosity in units of 10 42 erg/s. The 
result at a luminosity of L 42 = 100 with varying spectral index is shown for 
steep spectrum sources in Fig. 4 and in Fig. 5 for flat spectrum sources. The 
power law decrease can be observed up to a cut energy at E u ~ 10 9 GeV, 
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Fig. 4. Energy dependence of the generic 
AGN neutrino spectrum of steep spec- 
trum sources. The normalization is de- 
pendent on the particle index p, shown 
for p = 2.0, 2.2, 2.4, 2.6. The radio lu- 
minosity is assumed to be L42 = 100. 

where synchrotron losses of the neutr 
as the exponential cutoff of the spect: 



Fig. 5. Energy dependence of the generic 
AGN neutrino spectrum flat spectrum 
sources. The spectrum is shown at 
L42 = 100 and for three different parti- 
cle indices, p = 1.25, 1.50, , 1.75, 2.0. 

i-producing pion start to show as well 



4 Integration limits 



The integration over the luminosity has to be performed before the redshift 
integration, since the lower luminosity limit implicitly depends on the redshift. 
The limits for the z-integration are taken to be 



z m'm — 0.03 
^max 6 . 

The minimum of the z-integration is given by the fact that the flux is assumed 
to be diffuse. For z < 0.03, a contribution from the supergalactic plane is 
expected as will be discussed in section 5. The maximum redshift is taken 
to be z = 6, since the total z dependence decreases rapidly with redshift 
and any contribution above z max can be neglected as can be seen in Fig. 6 
(steep spectrum sources) and Fig. 7 (flat spectrum sources). Up to z = 6, the 
consistency of the use d models is ensured, since AGN have been detected up 
to redshifts of z « 6.4 (Willott et al.l O)]. 
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Fig. 6. Neutrino flux with variation of 
the upper redshift integration limit £ max 
for steep spectrum sources. The neu- 
trino energy is set to E® = 10 GeV 
while the particle index is varied as 
p = 2.0,2.2,2.4,2.6. 



Fig. 7. Neutrino flux with variation 
of the upper redshift integration limit 
Zmax for flat spectrum sources. The neu- 
trino energy is set to E® = 10 8 GeV 
while the particle index is varied as 
p = 1.25,1.50,1.75,2.00. 



The maximum radio luminosity limit is taken to be L^f x = 10 4 . Fig. 8 (steep 
spectrum sources) and 9 (flat spectrum sources) show the dependence of the 
flux at z = 2 and E° u = 10 8 GeV (flat) respectively E° = 10 GeV. The 
contribution is negligible for L.. l2 > 10 3 in the case of steep spectrum sources 
and for L 42 > 10 in the case of flat spectrum sources. This shows that most of 
the contribution to the flat spectrum neutrino flux comes from sources with 
moderate luminosities, L 42 < 10. The highest contribution to the flat spectrum 
neutrino flux comes from sources with luminosities around L 42 pa 300 — 1000. 



The lower radio luminosity limit depends on the following facts: 

(1) The absolute lower limit: For the model of extended sources, the jet-disk 
symbiosis is developed for FR-II galaxies and the absolute lower limit, 
L 42 min , is given by the lower luminosity limit for FR-II galaxies which 



is given as L i 78 pa 2.5 
approximately L 42 m i n = 



10 26 W/Hz with an integral radio luminosity of 
10. The sample of flat spectrum sources, on the 



other hand, ranges from 10 15 < L 42 < 10 4 and the absolute lower limit, 



is thus L{ 



10 



-1.5 



J 42,min 10 -^42,111111 

(2) The maximum proton energy which can be p roduced by an A GN is con- 
nected with the disk luminosity according to Lovelace! ( 197ft) : 



771 max 



C ■ y Lis 

up to a luminosity of L 46 pa 1, where the relation saturates, 



(8) 



jp max 



a. 



(9) 
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Fig. 8. Variation of the steep spectrum 
source neutrino flux with the upper in- 
tegration limit, using different spectral 
indices, p = 2.0, 2.2, 2.4 and 2.6. The 
redshift is arbitrarily taken to be z = 2, 
the energy is set to E® = 10 GeV. The 
flux normalization is saturated at ap- 



Fig. 9. Variation of the neutrino 
flux (flat spectrum sources), using 
p = 1.25, 1.5, 1.75 and 2.00. The 
redshift is chosen to be z = 2 and the 
energy is set to E® = 10 GeV. Sources 
at L > 10 do not contribute significantly 
to the flux. 



proximately L mSLX ~ 10 3 . 

The re sult is supported by the jet-disk symbiosis model, see Falcke and Biermarml 
( 19951 )]. Taking the highest observed energy so far, E p ps 10 21 eV and the 



maximum luminosity, L 



max 
disk 



10 erg/s, the constant C is determined 



to be C = 10 1L5 GeV. When the jet-disk symbiosis model by Falcke et al. 
is used and the ratio between E p and E u taken to be 20/1, the maximum 
neutrino energy is correlated to the radio luminosity due to the non linear 
relation between radio and disk luminosity, 



E] 
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This numerical relation is consistent with inserting all known parameters 
for the radio galaxy M87. C and (3 are given by the corresponding rela- 
tions between disk and radio luminosity for flat or steep spectrum sources 
and are listed in table 3. 

Therefore, the lower luminosity limit of an AGN contributing to a certain 
flux at a fixed neutrino energy is 



r mm 
-^42 



r j-s/f 



E 

I) >■ 



(10) 



The contribution to the flux for an energy corresponding to a lumi- 
nosity exceeding L disk = 10 46 erg/s, that is energies of E v = C/20 = 
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C [GeV] 


P 


El ut [GeV] 


steep 


2.2- 10 10 • (1 + z)- 1 / 4 


1/3 


~ 5 • 10 9 (> = 2) 


flat 


7.5 • 10 9 


0.3935 


~ 2 • 10 s 



Table 3 

Parameters for the determination of the lower luminosity integration limit. 
10 115 GeV^ 2 ■ 10 10 GeV is insignificantly small. 

Since the RLF is decreasing with luminosity, the energy spectrum will steepen 
from that energy on, when the luminosity-energy relation exceeds the value of 

the absolute luminosity limit, [-(f) > L^ min . This is the case for an energy 
E°*/GeV « {L s l f Mn f ■ C/GeV (11) 



in the comoving frame. The cut energies for the two models are given in table 3. 



5 Results 



5. 1 Discussion of the spectra varying the particle index 




10 12 

log(E/GeV) 




10 12 

log(E/GeV) 



Fig. 10. Neutrino spectrum of steep spec- 
trum AGN. The particle index is varied 
(p = 2.0, 2.2, 2.4, 2.6). 



Fig. 11. Neutrino spectrum of flat spec- 
trum AGN for different particle indices 
p= 1.25,1.50,1.75,2.00. 
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The final steep spectrum source flux is shown in Fig. 10. The spectrum is 
presented for four different particle indices, p = 2.0, 2.2, 2.4, 2.6. Here, p = 2 
is a commonly used value in literature and p = 2.6 has been computed in 
section 2 from the spectral index at 2.7 — 5 GHz for the steep spectrum popu- 
lation. The spectrum decreases as a power law up to a cut energy at approxi- 
mately 10 9 GeV. At higher energies, the spectrum is falling more rapidly due 
to pion synchrotron loss effects and the exponential cutoff of the spectrum. 
For a particle index of p = 2, the flux is approximately 

E 2 $ s « 10" 7 GeV/(s sr cm 2 ) 



for one neutrino flavor. This is very close to the current neutrino flux limit de- 
rived from the observed atmospheric muon and anti muon neutrino spectrum 



E 2 $ < 2.6 ■ HT 7 GeV/(s sr cm 



2> 



at energies around 300 TeV (Woschnagg and AMANDA Collaboration! (|2004f )]. 
The flat source neutrino spectrum is shown in Fig. 11 with a spectral index of 
p = 1.25, 1.50, 1.75 and 2.00. The cut energy lies at E ~ 10 9 GeV as discussed 
before. Values p ^ 2 are shown as a comparison, the final result will be given 
for p = 2 as indicated by diffuse shock acceleration. For p = 2, the flux is 



E 2 $ / « 10" 9 - 2 GeV/(s sr cm 2 ) 



for one neutrino flavor. This is about two orders of magnitude lower than the 
prediction for steep spectrum sources. 



5.2 Discussion of Uncertainties 



Apart from the lack of knowledge of the exact optical depth, the calculation of 
the neutrino flux $ bears uncertainties A$ because of two reasons: It is still 
a challenging task to make a precise prediction of the objects' redshift. Thus, 
the RLFs include high uncertainties for very distant objects. Nevertheless, 
the functions are quite well known up to a redshift of z ~ 2. Since objects 
at a higher redshift do not contribute significantly to the diffuse neutrino 
spectrum, the error due to uncertainties in distance measurements will be 
neglected compared to uncertainties in the jet-disk symbiosis model. 

The jet-disk symbiosis model agrees with the data when allowing a scattering 
of a factor of three in a L rai a (Ldisk) graph. This transfers to an uncertainty of 
A(p u = 0.78 ■ <j) u , that means the an uncertainty of a factor ~ 2 towards higher 
and a factor ~ 5 towards lower fluxes. The spectrum for flat spectrum sources 
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Fig. 12. AGN v^ + V^ neutrino spectrum 
with an index of p = 2. The flat 
spectrum sources (dotted line) do not 
contribute significantly compared to 
the steep spectrum population - the 
solid line represents the sum of flat 
and steep spectrum sources which does 
not differ from the result for steep 
spectrum sources. The data points 
result from unfolding the AMANDA 
neutrino spectrum. They follow the 
conventional atmospheric neutrino 
spectrum (dashed lines) and an upper 
limi t is derived as is indicated in the fig- 
ure I Woschnagg and AMANDA Collaboration! 



2 3 4 5 6 7 8 9 10 11 12 

log(E/GeV) 



Fig. 13. AGN v^+V^ neutrino spectrum 
with an index of p = 2.6 for steep 
spectrum sources and p = 2.0 for flat 
spectrum sources. The steep spectrum 
sources (dot-dashed line) make up 
most of the spectrum at lower energies 
(Eu > 10 6 GeV). The flat spectrum 
sources (dotted line) dominate the 
total spectrum (solid line) at energies of 
E u > 10 6 GeV. A signal that exceeds the 
atmospheric contribution is expected 
starting at E v w 10 8 GeV. 



includes an uncertainty of A(f>l = 0.9 • 
an a factor of 10 below the prediction. 



which means a factor of ~ 2 above 



Uncertainties due to the chosen cosmology would result in another factor of 
two in the change of the normalization. 

Since the maximum energy of the spectrum depends on the disk luminosity 
which is converted into a radio luminosity using the jet-disk symbiosis model, 
the cutoff may change by an order of magnitude at most. 



5.3 Discussion of a source flux from the supergalactic plane 



It has been stated before that the calculations include sources at z > 0.03 
to exclude possible anisotropics from the supergalactic plane. In the following 
paragraph, a possible flux from the supergalactic plane will be discussed. Note 
that this is just to make an estimate of a possible contribution from nearby 
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sources and that final results do not include sources at z < 0.03. Assuming 
that the contribution to the total flux at z < 0.03 is concentrated in the 
supergalactic plane yields an inc r ease o f the flux by a factor 10 since the 



supergalactic plane, Lahav et al. ( 2000f )]. covers approximately 10 % of the 



sky. The ratio of the diffuse steep spectrum flux at z > 0.03 and the flux from 
the supergalactic plane is of the order 



$steep( z < o 03) 
$ stee f(> > 0.03) 



0.01 



Since this is only some percent of the total diffuse flux, it is negligible. For flat 
spectrum sources, the result is 

v L ~ n 02 

$flat( z > Q.03) 



and thus also negligible compared to the flux from outside the supergalactic 
plane. 



5.4 Discussion of the final spectra 



The comparison of the two source type spectra considering only one neu- 
trino flavor 1 in Fig. 12 shows that the dominant contribution comes from 
steep spectrum sources assuming a particle index of p = 2.0. However, us- 
ing p = 2.6 for the steep source spectrum reduces the contribution from this 
sample significantly as is presented in Fig. 13. The cutoff in the energy spec- 
trum at E v ~ 10 9 GeV occurs since the maximum e nergy for each s ingle 
source depends on the sour ce luminosity as shown by Lovelace! ( 19761 ) and 
Falcke and Biermann ( 199.4 ) . Thus, the cutoff energy of the isotropic dis- 



tribution depends on the sources' luminosity distribution. The atmospheric 
neutrino spectrum (dashed lines) dominates the spectrum up to energies of 
~ 10 7 GeV. 



The maximum optical depth of the sources can be estimated by assuming that 
the contribution of the calculated spectrum has to be smaller than the limit 
set by AMANDA mentioned above. Thus, in case of an E~ 2 spectrum for steep 
spectrum sources, the optical depth can be restricted to r e // < 2. This factor 
will grow significantly when using an E~ 2S spectrum, i.e. r e // < 1000. For flat 
spectrum sources, an effective opacity of r e // < 400 can be estimated. 

1 The total flux is divided by a factor of three, taking into account neutrino oscil- 
lations. 
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A comparison of the results using two different spectral indices for the steep 
spectrum sources makes clear how much the total flux contribution depends 
on the spectral index of the sources. Changing the spectral index from p = 2.0 
to p — 2.6 reduces the contribution from the source sample by two orders of 
magnitude. In the pessimistic case, the detection of a signal from these types 
of AGN will require an array like IceCube or Auger which are sensitive to 
UHE neutrinos at ~ 10 19 eV. In the case of an E~ 2 spectrum the detection of 
a diffuse signal should be possible in the near future. Fig. 14 shows a summary 
of var ious extragalactic neutrino flux predictions, see (Learned and Mannheim 
( 2000| )] for a review. In the usual E~ 2 scenario, the AGN spectrum should be 
detectable in the very near future, presumably already by running neutrino 
experiments like AMANDA. If the spectrum is steeper, however, it has been 
shown that a detection of a significant neutrino signal from AGN would need 
some years of operation of a km 3 experiment like IceCube. 
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Fig. 14. + z7 M neutrino flux predictions. (1) and (2) are the predic- 
tions done in this paper, BBR(I) using p = 2.6 for the steep spectrum 
source s, BBR(II ) applying y = 2; (3): pj interactions, MPR Mannheim et al.l 
1 2001)]: (4): M [Mannheiml lll995| )]. pj in blaz ar and yy in its host galaxy ; (5): 
RB |Rachen and Biermamll|l j?7; (6): WB |Wa,xma,Ti and Bahcalll |l99flh ]. dif- 
fuse neutrino flux from GRBs, assuming that GRBs are the sources of cosmic rays 
at E > W 19 eV. The dashed lines indi cate the atmospheric neutrino spectrum ac- 
cording to Volkova and ZatsepirJ l|l980l) . the upper line representing the horizontal 
flux, the lower line showing the vertical contribution. The do tted lines rep r esent a 
model for the atmospheric charm contribution as predicted in (Martin et al. 1 2003 )]. 
The uncertainties in the model occur due to larger errors in measurements of the 
Bjorken variable in forward scattering. Data points are AMANDA data showing the 
unfolded spectrum of 2000 data. It follows the atmospheric prediction. The limit on 
the extraterrestrial neutrino flux is set based on the data from the last bin of the 
unfolded spectrum Woschnagg and AMANDA Collaboration! (|2004h ]. 



21 



6 Acknowledgments 



The authors would like to thank K. Mannheim and H. Falcke for valuable 
discussions which led to much progress in realizing this work. 

The work with P. L. Biermann has been supported through the AUGER theory 
and membership grant 05 CU1ERA/3 via DESY/BMBF (Germany). Further 
support for the work with P. L. Biermann has come from the DFG, DAAD, 
Humboldt Foundation (all Germany), grant 2000/06695-0 from FAPESP (Brasil) 
with G. Medina- Tanco, a grant from KOSEF (Korea) through H. Kang and 
D. Ryu, a grant from ARC (Australia) through R. J. Protheroe, and a Euro- 
pean INTAS grant with V. Berezinsky. Special support comes from the Eu- 
ropean Sokrates/Erasmus grants in collaboration with East-European Uni- 
versities, with partners M. Ostrowski, K. Petrovay, L. Gergely, A. Petrusel, 
and M.V. Rusu, and VIHKOS through the FZ Karlsruhe. Current support 
comes from NATO for a collaboration with S. Moiseenko and G. Bisnovatyi- 
Kogan (Moscow), and the Chinese Academy of Sciences of China in work with 
Y. Wang (Beijing/Nanjing). 



22 



References 



Biermann, P. L., 2001. High energy particles in Active Galactic Nuclei. In: Shapiro, 
M. M., et al. (Eds.), Kluwer Acad. Publ. p. 115, review at the Erice school Nov. 
2000. 

Biermann, P. L., Strittmatter, P. A., 1987. apj 322, 643. 

Bridle, A. H., Perley, R. A., 1984. Ann. Rev. Astron. Astrophys. 22, 319. 

Dunlop, J. S., Peacock, J. A., 1990. Month. Not. Roy. Astr. Soc. 247, 19. 

Falcke, H., Biermann, P. L., 1995. Astron. & Astrophys. 293, 665. 

Falcke, H., Malkan, M. A., Biermann, P. L., 1995. Astron. k Astrophys. 298, 375. 

Gregorini, L., et al., 1984. Astron. J. 89, 323. 

Hillas, A. M., 1984. The Origin of Ultra-High-Energy Cosmic Rays. Ann. Rev. As- 
tron. Astrophys. 22, 425. 
IceCube Collaboration, 2004. Astropart. Physics 20, 507. 
Kiihr, H., et al., 1981. Astron. J. 86, 854. 
Lacy, M., et al., 1999. Month. Not. Roy. Astr. Soc. 308, 1096. 
Lahav, O., et al., 2000. Month. Not. Roy. Astr. Soc. 312, 166. 

Laing, R. A., Riley, J. M., Longair, M. S., 1983. Month. Not. Roy. Astr. Soc. 204, 
151. 

Learned, J. G., Mannheim, K., 2000. Annu. Rev. Nucl. Part. Sci. 50, 679. 
Lovelace, R. V., 1976. Nature 262, 649. 
Mannheim, K., 1995. Astropart. Physics 3, 295. 

Mannheim, K., Protheroe, R. J., Rachen, J. P., 2001. Phys. Rev. D 63, 23003. 
Martin, A. D., Ryskin, M., Stasto, A., 2003. Acta Phys. Polon. B34, 3273. 
McGilchrist, M. M., Riley, J. M., 1990. Month. Not. Roy. Astr. Soc. 246, 123. 
Miyaji, T., Hasinger, G., Schmidt, M., 2000. Astron. & Astrophys. 353, 25. 
Montaruli, T., 2003. In: Particle Astrophysics Instrumentation. Edited by Peter W. 

Gorham. Proceedings of the SPIE, Volume 4858. p. 92. 
Miicke, A., et al., 1999. Publ. Astron. Soc. Aust. 16, 160. 
Miicke, A., et al., 2000. Comput. Phys. Commun. 124, 290. 
Particle Data Group, 2004. Phys. Let. B 592, 1, Particle Physics Booklet. 

URL |http : //pdg . lbl . gov| 
Pauliny-Toth, I. I. K., Kellermann, K. I., 1972. Astron. J. 77, 797. 
Pauliny-Toth, I. I. K., et al., 1972. Astron. J. 77, 265. 
Pauliny-Toth, I. I. K., et al., 1978. Astron. J. 83, 451. 
Peacock, J. A., 1985. Month. Not. Roy. Astr. Soc. 217, 601. 
Peacock, J. A., Wall, J. V., 1981. Month. Not. Roy. Astr. Soc. 194, 331. 
Protheroe, R. J., 1998. In: Towards the Millennium in Astrophysics, Problems and 

Prospects. International School of Cosmic Ray Astrophysics 10th Course, p. 3. 
Rachen, J. P., Biermann, P. L., 1993. Astron. & Astrophys. 272, 161. 
Rawlings, S., Eales, S., Lacy, M., 2001. Month. Not. Roy. Astr. Soc. 322, 523. 
Rybicki, G. B., Lightman, A. P., 1979. Radiative processes in astrophysics. J. Wiley 

& Sons, Inc. 
Schechter, P., 1976. Astrophys. J. 203, 297. 
Spergel, D. N., et al., 2003. Astrophys. J. Sup. S. 148, 175. 
Volkova, L. V., Zatsepin, G. T., 1980. Sov. J. of Nuc. Phys. 37, 212. 
Wall, J. V., Peacock, J. A., 1985. Month. Not. Roy. Astr. Soc. 216, 173. 



23 



Waxman, E., Bahcall, J., 1999. Phys. Rev. D 59, 23002. 

Willott, C. J., McLure, R. J., Jarvis, M. J., 2003. Astrophys. J. Let. 587, L15. 
Willott, C. J., et al, 2001. Month. Not. Roy. Astr. Soc. 322, 536. 
Woschnagg, K., AMANDA Collaboration, 2004. Astro-ph/0409423, talk at Neutrino 
2004. 

Wright, A. E., et al., 1982. Austr. J. of Phys. 35, 177. 



24 



